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Using the simplest possible ingredients of a rupture model with thermal fluctuations, we provide 
an analytical theory of three ubiquitous empirical observations obtained in creep (constant applied 
stress) experiments: the initial Andrade-like and Omori-like 1/t decay of the rate of deformation 
and of fiber ruptures and the l/(t c — t) critical time-to-failure behavior of acoustic emissions just 
prior to the macroscopic rupture. The lifetime of the material is controlled by a thermally activated 
Arrhenius nucleation process, describing the cross-over between these two regimes. Our results give 
further credit to the idea proposed by Ciliberto et al. that the tiny thermal fluctuations may actually 
play an essential role in macroscopic deformation and rupture processes at room temperature. We 
discover a new re-entrant effect of the lifetime as a function of quenched disorder amplitude. 
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Constant stress (so-called "creep") experiments con- 
stitute a standard testing procedure in material sciences. 
The typical response to the sudden application of a con- 
stant stress is that the strain rate as well as the acoustic 
emission rate first jump to high values followed by slow 
universal power law decays, respectively called the An- 
drade law 0] for the strain rate and the Omori law Q 
for the acoustic rate. Then, after a long decay whose du- 
ration may vary within extraordinary large bounds (see 
below), the rates rebound and accelerate (while the ap- 
plied stress remains constant) by following a power law 
acceleration resulting in a finite-time singularity (the rup- 
ture of the sample) . The two regimes of decelerating fol- 
lowed by accelerating rates and the lifetime of the struc- 
ture are the result of a subtle interplay between the pre- 
existing micro-heterogeneity of the material and the self- 
organized evolving deformation and damage due to dislo- 
cation motion and/or micro-cracking. Up to now, there 
are no theory encompassing all these regimes. Here, we 
propose a simple mechanism that provides an explana- 
tion of all these observations, which is based on the recent 
proposal 

0,0 

that thermal noise is strongly renormal- 
ized by quenched heterogeneities. Based on the analy- 
sis of a simple fiber bundle rupture model, Refs.0, § 
showed that the average lifetime of the fiber-bundle takes 
an Arrhenius form with an effective temperature renor- 
malized from the bare temperature T to a value strongly 
amplified by the presence of the frozen disorder in the 
rupture thresholds f c (i), hi agreement with experiments 
and numerical simulations. This result suggests that the 
usual assumption of neglecting the role of thermal fluctu- 
ations in material rupture processes at room temperature 
may actually be incorrect (see [f| for early discussions): 
due to frozen heterogeneities, tiny thermal fluctuations 
can be amplified many times, thus actually controlling 



the time-dependent aspects of failure. Our purpose is 
to extend the analysis of this model by showing that it 
is able to reproduce all the empirical observations men- 
tioned above in creep (constant applied stress) experi- 
ments. This shows that the simplest possible ingredients 
of a rupture model together with thermal fluctuations can 
render essentially all of the richness of creep experiments. 

The democratic fiber-bundle model (DFBM) with 
thermal noise 0, 0] can be seen as a mean field treat- 
ment of rupture. A macroscopic constant load F — Nf 
is applied at time t = to an initially undamaged system 
made of a very large number N of parallel elastic fibers 
(the results derived below are obtained in the thermody- 
namic limit N — > oo). At all times, F is shared democrat- 
ically among all (1 — $(t))AT surviving fibers, where Q(t) 
is the fraction of broken fibers at time t. The externally 
applied force per surviving fiber is thus 



fa 



fa 



1 - $(t) 
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The strength of each fiber i is characterized by a criti- 
cal value f c (i) drawn for a distribution Pd{f), centered 
on the mean equal to 1 and with variance T<j. Putting 
the mean strength to 1 sets the force scale. The hetero- 
geneous strengths are constant and determined once for 
all, corresponding to a frozen disorder, which is "read" 
in a certain organized way as the rupture develops. Mi- 
croscopic thermal fluctuations are taken into account 
by assuming that a fiber with load f a and threshold 
fc(i) > fa has a non-zero probability G(f c (i) — f a ) to 
rupture per unit time governed by the rate with which 
a thermal fluctuation can activate a microscopic force 
A/j > f c (i) — f a to pass the rupture threshold f c (i): 

G{f c (i) - f a ) = 2 erfc ( ^v^^ )' wn ere erfc(x) is the 
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complementary error function, T is the variance of the 
thermal force fluctuations A/, and 7 is a microscopic con- 
stant rate fixing the time scale of the thermal activation 
process. This expression amounts to introducing a zero- 
mean normal distribution of thermal fluctuation forces 
A/j with variance T and with correlation time propor- 
tional to I/7. 

We first follow |j] and introduce the distribution 
Q(f, t) of the rupture thresholds of the unbroken fibers at 
time t. Obviously, Q{f,t) — for / < f a , since all these 
fibers are already broken. Politi et al. Q| have shown 
that Q(f, t) can be approximated with a very high accu- 
racy in the limit ./V — > 00 by the initial distribution Pd(f) 
of rupture strengths truncated at a lower value f s (t) , 



Q(f,t) = Pd(f) , for / > f„(t) 



(2) 



and otherwise, where f s (t) is determined by the self- 
consistent equation 



r/ s (t) 

$(*) = / df Pd{f) 



(3) 



expressing that all fibers whose strengths are below f s (t) 
have failed at some time before t. This approximation for 
Q(f, t) with <|3J) amounts to view the time-dependent rup- 
ture as a "front" propagating and "eating" the distribu- 
tion Pd(f) from the weakest towards the strongest fibers. 
We also have by definition <I>(t) = 1 - df Q(f,t). 
Taking the time derivative of <E>(i) and replacing Q(f,t) 
by —Q(f, t) G(f — f a ) expressing that the rate of break- 
ing is controlled the thermally activated rupture pro- 
cess acting on each fiber independently, we get $ — 
J-^ d f Q(f,t) G(f ~ fa)- Putting © in this equation 
and taking for Pd{f) a normal distribution centered on 1 
with variance Td as in 0, yields 



$ = 



2Jf. y/2^T d 
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Making explicit Pd(f) in 10 gives 



df . 

(4) 



$ = 



erf 



/.-I 



, /, = 1 + ypZfi irf (2$ - 1), 

(5) 

where y = irf (z) is the inverse function to the error func- 
tion z = erf (y). Putting f s in J3J| gives 



<E> = R(4>) = ±J^ erfc [£(*, z)] dz , (6) 



L($,z) = 



f - 



1 - $ 



+ /xirf(22f-l), (7) 



with = ^/TdjT. The solution of equation ijfjj with Q 
provides in principle all the information on the fraction 



$(t) of broken fibers. This equation (JTJ is valid as long 
as the approximation holds (see below). 

The parameter fi quantifies the relative importance of 
the thermal fluctuations compared with the quenched 
heterogeneities. The relevant regime for applications to 
macroscopic ruptures at room temperature is /i > 1 and 
often > 1, that is, thermal fluctuations are tiny con- 
tributions to the applied macroscopic mechanical forces. 
Indeed, assuming that the energy barrier to rupture a 
fiber corresponds to the Griffith energy ss g c 2 neces- 
sary for nucleating a crack of half-length c in the solid 
with surface energy g, we obtain /i sa 1.5 — 4 • 10 3 for 
c = 1 micron and /j, ~ 1.5 — 4 for c = 1 nanometer, using 
g = 10 — 50 erg/cm 2 for most solids. Thus, even for the 
smallest microcracks, thermal fluctuations are very small 
in relative value. 

It turns out that this regime fi > 1 allows for a very 
convenient approximation of -/?,($) obtained by lineariz- 
ing z) with respect to z. Then, the integral over z 
in JHJ) can be calculated explicitly to yield 



$ = i?($) 
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where 



u($) = ri 2 ($,$) = - [/,($) - f a (<s>)} 2 



(8) 



(9) 



and £>($) = {l/v^T d )df s {<S>)/d<5> = exp{irf 2 (2$ - 1)}. 
Notice that £/($) in © has a clear physical interpreta- 
tion. It is the energy barrier between the actual force 
/a($) O and the front value / s ($) © of the distri- 
bution Q(f,t) in @. Equation JSJ) is valid as long as 
U($) » T, which implies $ < <f> c , where $ c is such 
that the force f a per surviving fiber reaches the aver- 
age strength 1: / a ($ c ) = 1 yielding $ c = 1 - f . Such 
fractions $ — > $ c correspond to the ultimate regime of 
explosive failure. We have checked by direct numerical 
calculations that JSJ) provides an exceedingly precise ap- 
proximation of equation (jBJ) as long as /x > 1 and not too 
close to $ c (in practice to within a few percent). 

To go further, we need to distinguish between two 
regimes, <£><$* and $ > where <&*, solution of 
dR(&)/d$ = 0, corresponds to the minimum failure rate. 
For fj, > 1, $* is actually independent of the temperature 
T and is the root of the equation 



£>($*)(! - $* 



fo 



V^TTTd 
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where a is a important physical parameter quantifying 
the strength of the disorder relative to the applied force. 
The explicit approximate solution of equation 1101) is 



**(a) 




a < 3/2, 
a > 3/2. 



(11) 
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For example, this gives <&*(a = 2) = 0.089 compared 
with the exact value 0.092. 

It follows from (JSJ that the time to reach some <& 
is given by -yTt = 4tth /* D(z)U(z) e u ^' T dz. For 
< $ < $*(a), due to the exponential factor e u / T , 
the main contribution to the last integral comes from a 
small neighborhood of the upper integration limit. This 
yields 



7 t ~ 47T/U 



£>($)[/($) 



,!7(*)/T 



(12) 



for $<<!>*. This approximation is correct under the 
assumption that e u ^^ T is rapidly increasing with $, 
i.e., if |j4($)| > T. The absolute value |...| stresses 
that this condition applies also for «3> > $*(a). The same 
reasoning in this case gives a similar approximation 
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where £ c — i is the time to complete rupture. The con- 
dition \A(&)\ > T shows that both relations and 
(I13JI do not work in the vicinity of the minimum rate 
of fiber failures given by the solution of 1)11) |). for which 
A($*) = 0. For < $ < $*, combining © and JUJ, 
we obtain U = T/t for t < t* , where t* is defined by 
This gives 



[/[$(*)]= Tin 7 t (t<i*) 



(14) 



valid for 7i ^> 1 (this condition simply means that the 
thermal fluctuations have had time to contribute several 
independent jolts). The constant of integration gives the 
In 7 contribution determined from matching with the ini- 
tial stage. Replacing the lhs of © by (|14|) and putting 
f a ~ fo (for <E> small) gives, in view of equation (JSJ), the 
fraction rate 



47r/ii hi^ft 



exp 



1 

2^ 



1 - f - v / 2rin 7 t 



(15) 

Expression 1)15(1 is one of our main results: the failure 
rate $ of fibers decreases after application of the load 
proportionally to l/t, up to logarithm corrections. This 
l/t decay lasts as long as $ remains smaller than $*. 
This l/t law is known in seismology as the Omori law. 
It is also ubiquitous in creep experiments with expo- 
nents that are often close to or smaller than our pre- 
diction 1. For intermediate times such that "ft < e 1 / 21 " , 
~ t S~l e( 1 ~f°^ 2T ^Tt/Td^ w hich gives an apparent 
exponent - 1/& with p < 1. For lnT* > (1 - f ) 2 /2T, 
p — ► 1 + (T/Td) which is close to but slightly larger than 
1. Numerical simulations confirm these predictions ac- 
curately. See for instance figure 2 of which our the- 
ory explains quantitatively. Exact numerical integration 
and our analytical approximation coincide everywhere, 
excluding a time interval corresponding to a very small 



vicinity of the stationary point $*. Note that Andrade's 
law also derives from the deformation rate being pro- 
portional to df a (t)/dt = / $/(l - $(i)) 2 oc $ as 
varies much more slowly than $. 

Let us now turn to the description of the second regime 
> $*, relevant to obtain the failure rate up to global 
failure. Combining (JSJ and l|13|) . we obtain the expression 



U [$(*)] = T ln[7(* c - t)] ($*<$< $ c 



(16) 



which is analogous to 1(11)1. The regime $* < < <1> C 
is strongly influenced by thermal fluctuations, so that the 
disorder term can be neglected to obtain, in view of i|16|) 

and ©, $ c 



$(*) 



/o^piZSHL. Differentiating 

JU \~sJlT ln 7 (t c -t) b 

both sides of this expression with respect to t yields the 
failure rate 



*(t) = C(t)/(t c - t) , 



(17) 



where C(t) = f T/[c(l - c) 2 ] with c = ^2T\n[-y(t c - t)}. 
This is the second important result of our analysis (see 
also eq.(Bll) in 3]), which shows that, for $ > $*, the 
failure rate accelerates towards a finite-time singularity 
approximately as ~ l/(t c — t). Such a behavior has been 
documented extensively in experiments on rupture of het- 
erogeneous material 6J. Our analysis provides a novel 
mechanism for the ubiquitous time-to-failure regime ob- 
served in heterogeneous material. Strong quenched het- 
erogeneity has been shown to play an essential role in 
controlling the critical nature of the rupture process Q 
and in the existence of a time-to-failure power law such 
as <|17[) . Here, we confirm that the heterogeneity is essen- 
tial to renormalize the thermal fluctuations While 
the philosophy is similar, the mechanism is different and 
novel. As for the Omori law, the logarithmic corrections 
in (|17[) may give an apparent exponent of the power law, 
slightly smaller than 1, as observed in experiments. Our 
numerical tests show that expression (|17|l provides an ap- 
proximation which coincide almost everywhere with the 
exact solution inside the interval t* < t < t c . 

There is a simple physical interpretation of the tran- 
sition between the two above mentioned rate behav- 
iors 1)15)1 and l|17jl. To explain the first (rate decaying) 
regime, consider the degenerate case of spontaneous frac- 
ture (f = 0, $* = 1) for which [/($) = / 2 /2. As time 
increases, f s grows, the remaining fibers are stronger and 
the failure rate decays together with the rate of change 
of the energy barrier. The second regime can be quali- 
tatively understood by taking the limit of zero disorder 
(T d = 0, $* = 0), leading to U(§) = (1 - / a ) 2 /2. The 
force f a per remaining fiber grows with time, the fibers 
break more and more easily and the failure rate grows to 
give the fracture in finite-time. In the intermediate case 
< <f>* < $ c , due to the competition between the growth 
of f s and / a , the two regimes co-exist. At early times, the 
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growth of f s dominates giving the Omori and Andrade 
laws, followed by the growth of f a in the second regime 
$ > $* giving the power law finite-time singularity. 

Lastly, we turn to the behavior for $ « <£>*, which 
turns out to provide the dominant contribution for the 
total time for rupture, as shown in 0, U Indeed, the 
fiber bundle spends most of its time in the vicinity of 
the stationary point <&*, corresponding to the minimum 
failure rate. In this case, U can be expanded as 



t/($) = [/($*) -£(<!>*)($ -$*) 2 , £($) = - 



and equation © becomes 



T 



$ ~ i?(<f>*) exp 
The solution of this equation is 

erfi 



ld 2 [/($) 
2 d$ 2 : 



(18) 
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= 2 



B($*) 

7TT 



#($*)(*-**), 
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where erfi (z) = 4 erf (iz) is the imaginary error function. 
Using its asymptotics erfi (z) ~ -^=^ e z for large z to- 
gether with 119|) . equation (|18l) becomes 



T 



rf* 2 

~dT ~ B($*) (t-t* 



\[f = $ - $* 



whose solution yields the fracture rate 
V-sign (t-i*) 



2|t-t*|- v /B($*) ln(7|t-t* 



(7|*-t*| »1) • 



Expression (|18|l allows us additionally to calculate the 
total lifetime t r of the fiber bundle: 



7*c 



1 



exp 



B($*) , ^2 

^r 1 (<J>-$*) 2 



The calculation of this integral with the use of 

[/($* 



(20) 



7 t c ~ 4tt0F 7== ~frT~ exp 



T 



(21) 



This expression, together with 1[11[) . recovers the main 
result of 0, 0] , while improving on the prefactors to the 
main Arrhenius-type dependence. 

Using © and ((TUjl . [/($*) can be written explicitly 




(22) 



where the sign + (resp. — ) corresponds to the case 
$* > 1/2 (resp. $* < 1/2). As shown in Fig.HJ E/($*) is 
a non-monotone function of T^. Due to the above men- 
tioned competition between quenched disorder and the 



growth of the actual force f a , [/($*) decreases as long 
as Td < and then increases with increasing Td be- 
yond . The first regime Td<T d * corresponds to the 
effect discovered in Refs. 0, @ and mentioned above of 
the renormalization of thermal fluctuations by quenched 
disorder, and consequently of decreasing strength by in- 
creasing the disorder. Since a larger £/($*) corresponds 
to a large lifetime through 1)21(1. we uncover the new ef- 
fect of a strengthening of the fiber system by increasing 
the disorder beyond a certain threshold. All our formulas 
have been checked by direct numerical integration with 
excellent agreements. We expect that extensions of the 
DFBM to non-mean field power law interactions [8j will 
not change our results qualitatively but may modify the 
Omori's and time-to-failure exponents. 
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FIG. 1: Effective barrier energy [/(<&*) as a function of the 
disorder strength Td, in the case where $* = 1/2, correspond- 
ing to T d * = | / 2 , for f = 0.1 leading to T d * ~ 0.025. This 
illustrates the non-monotonous behavior of [/(<!>*) and thus 
of the lifetime t c with Td. 
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